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Abstract. As a population, field massive stars are relatively enigmatic, and this re- 
view attempts to illuminate this sector of the high-mass stellar population, which com- 
prises 20 - 25% of the massive stars in star-forming galaxies. The statistical properties 
of the field population are vital diagnostics of star formation theory, cluster dynamical 
evolution, and stellar evolution. We present evidence that field massive stars originate 
both in situ and as runaways from clusters, based on the clustering law, IMF, rotation 
velocities, and individual observed in situ candidate field stars. We compare the known 
properties of field and cluster massive stars from studies in the Magellanic Clouds and 
the Galaxy, including our RI0TS4 complete spectroscopic survey of SMC OB stars. In 
addition to the origin of the field massive stars, we discuss additional properties includ- 
ing binarity, runaway mechanisms, and some evolved spectral types. 



1. Introduction: Definition and frequency of field OB stars 

In the last session, we discussed massive stars within clusters, and I now turn to a more 
enigmatic component of the high-mass stellar population: the field stars. What is the 
nature of the field massive star population? In particular, are there OB stars that form 
in situ in the field, more or less isolated? Or are most massive field stars runaways that 
are ejected from clusters? Note that "field" stars are subject to different definitions, es- 
pecially regarding the degree of isolation, both in the number of immediate companion 
stars and distance from massive clusters. Thus what follows is a review of studies that 
adopt necessarily heterogeneous definitions of the field. 

For a simple definition of field OB stars, we consider the clustering law, which is 
the number of OB stars A^* per cluster. This is simply a manifestation of the cluster mass 
function, and both are known to robustly follow a -2 power-law distribution, which for 
the clustering law is: 

A^(A^*) dN, oc dN, . (1) 

If we consider only the OB stars and ignore all the lower-mass stars, then we find that 
this power-law extends s moothly down to A^* = 1, i.e., individual OB stars, which 
we define to be field stars dOey et aLll2004 i). We obtained this result t ising the UBVR 



photometric survey of the entire Small Magellanic Cloud (SMC) by iMasseyl ([2002), 
uniformly identifying all OB star candidates from the criteria that the reddening-free 
parameter Qubr < -0.84 and B < 15.21. Following fBattinellil (Il99ll) . we apphed a 
friends-of-friends algorithm to assign stars to clusters, adopting the clustering length 
that maximizes the number of identified clusters, which is 28 pc. OB candidates having 
no other OB candidate within this distance are then defined to be field stars. 
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That these field stars are simply a smooth extension of the power-law clustering 
suggests that they are dominated by objects that are simply the "tip of the iceberg" 
on small, sparse clusters containing only the single OB stars that formed in situ as 
members of these groups. However, note that since the clustering law is presented as 
a logarithmic relation, we can also accommodate a significant contribution, even 50%, 
from runaway stars, and still present a smooth power law for the clustering law. 

From this power law (equation [B, the fraction of field massive stars is (Oey et al. 
2004): 

/(field) - (lnA^*,up + 0.5772)"' , (2) 

where A^«,up is the number of OB stars in the richest cluster in the entire galaxy or 
ensemble. For most astrophysical situations, /(field) is in the range of about 1/3 to 
1/10. For the SMC, equation |2] predicts /(field) ~ 0.20, while the observed value is 
about 0.26. This is com parable to the estimated frequency of ~20% for O stars in the 
Milky Way (iMason et al...2009. : .Gies.l987,) . 



2. The Field IMF 

Given this uniformly identified sample of field OB candidates for the entire SMC, we 
are in the final stages of a survey to obtain spectroscopic observations for all of these 
stars: the Runaways and Isolated 0-Type Star Spectroscopic Survey of the SMC (RI- 
0TS4). This is being carried out at the Magellan Baade 6.5-m telescope by multislit 
observations with the IMACS imaging spectrograph at a resolution of R ~ 2600. This 
yields spectroscopic classifications and radial velocities good to 10 - 15 km s"'. Our 
coverage of the SMC bar region is essentially complete, while the SMC wing region is 
currently about 85% complete. 

A primary goal of our survey is thus to provide a definitive measurement of the up- 
per IMF in the field, which is presented a t this meeting by Joel Lamb. As is well known. 



Massey (2002) and iMassey et all (119951) found an upper IMF power-law slope that is 
much steeper than the Salpeter value F = -1.35 in both the Large and Small Magel- 
lanic Cloud from observations of limited areas in these galaxies, obtaining F 4. Our 

preliminary results confirm a steeper slope, though with a less ex t reme value: F 3 



(Lamb et al. 2011, these proceedings). Similarly. Ivan den Berghl (120041) evaluated the 
distribution of Milky Way O stars as a function of spectral type and found a preponder- 
ance toward late types in the field relative to associations. This is thus consistent with 
a steeper field slope in our Galaxy as well. 

Since Wolf-Rayet (WR) stars have higher-mass progenitors than typical OB stars, 
we can also examine the relative frequency of WR stars to probe the field and cluster 
IMF. In the SMC, t wo WR stars are in our field sample, out of a total of 12 known WR 
stars in this galaxy (Massev & Duffv"200ll: lMassey et atl bOOlV This yields a field WR 



frequency of 17% relative to 26% for OB stars. While this difference is not statistically 
significant because of the tiny number of WR stars, we note that the results do suggest 
a lower WR frequency, as expected for a steeper field IMF. 



3. Oe, Be, and B[e] Stars 



Classical Oe and Be stars are another type of massive star that are easy to identify from 
their strong Balmer emission lines. However, the origin and frequency of this phe- 
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nomenon is less straightforward than for WR stars, since apparently only some fraction 
of OB stars go through an Oe/Be phase. Field vs cluster properties for these stars are 
best examined in the Magellanic Clouds, but no clear differences have been identified in 
the distributions of spectral type, luminosity class, and rotation velocities dKeller et al. 



19991 : ^Martava n et"al]|2006l l2007h . Evidence for variation in their frequency between 



field and clusters is also inconclusive. Our RI0TS4 survey yields a preliminary fre- 
quency of 53% and 21% for field Be and Oe stars, respectively. These high frequen- 
cies reflect the known prevalence of Oe/Be stars at low metallicity for the SMC (e.g., 
iMaeder et al.. 1999 : Martayan et al.. 2006 ). 

We identi fied a new B [ el supergiant in the course of the RI0TS4 survey, the star 
LHA115-S62 (iHenizel 119561) . Thus we have one B[e] supergiant in the field relative 
to the total of 6 in the SMC, adding to the five previously known B[e] supergiants 
(IWisniewski et al.ll2007h . This yields a frequency of 17% B[e] supergiants in the field, 
which is the same value found for WR stars. 



4. Binary Stars 

Statistics for binary stars are an especially important diagnostic for the origin of field 
stars, since theories for both star formation and dynamical ejection from clusters can 
make specific predictions for the frequencies and binary parameters for field stars. Be- 
cause of the difficult and time-consuming nature of binary identification across the 
vast parameter space in binary mass ratio, separation, inclination, and period, differ- 
ent techniqu es are needed to probe the full parameter space (e.g., Sana & Evans 201 H; 
iMason et aLll2009l) . ranging from high-resolution imaging with speckle interferometry 
to high-resolution spectroscopic monitoring. Furthermore, obtaining reliable statistics 
also depends on observing complete samples over appropriate time domains. Mason 
et al. (2009) have compiled the most complete statistics for massive star binaries to 
date, examining 347 Milky Way O stars in a magnitude-limited sample catalogued by 
|Maiz-Apellaniz et al. (2004). To date, they find that the field and cluster O stars have 
frequencies of 59% and 75%, respectively, that are at least binary or with greater mul- 
tiplicity. While this difference is statistically significant, it is important to bear in mind 
that binary frequencies for both populations are lower limits. Sana & Evans (201 1) find 
a cluster binary frequency around 50% for 9 Galactic OB associations; since these are 
typically farther than the stars in the Maiz-Apellaniz et al. sample, the binary detection 
rate is lower. The binary frequencies for thei r extragalactic clust ers are lower still. The 
frequency of binary runaway O stars is 43% dMason et al.ll2009l) . 

Bragan9a et al. (2011, in preparation) obtained observations of 379 solar neigh- 
borhood B stars with the MIKE echelle spectrograph at Magellan. Preliminary results 
show that about 21% of these stars are double-lined spectroscopic binaries, with anal- 
ysis of field vs cluster members in progress. Our RI0TS4 survey of the SMC includes 
binary monitoring of three fields, corresponding to about 10% of the sample. With ob- 
servations of a dozen epochs in hand, about 9 out of 16 non-Be stars are single-lined 
spectroscopic binaries, corresponding to 60% of the sample. We caution that stochas- 
ticity is an important factor in this preliminary result. 

As mentioned above, the statistics for binary stars promise to be critical diagnos- 
tics for the origin of both in situ and runaway field stars. At present, reliable statistics for 
binary orbital parameters have not yet been published, but we stress the importance of 
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obtaining these data. For example, we almost have useful estimates of the non-compact 
binary runaway frequency, which is a critical probe of dynamical ejection (§6). 



5. Rotation Velocities 



Rotation velocities are the parameter for which the field and cluster stars show the most 
revealing differentiation. It has long been known that field stars show lower values 
of vsin/ than cluster stars (e.g., .Guthrie .1984) . This difference is a critical clue to 
the origin of field stars, and the proposed models reflect the different possible mecha- 
nisms for these origins. IWolff et aH (l2007l l2008h suggest that the different rotational 
velocities result from the different star-formation environments, with lower turbulence 
in small star-forming regions resulting in lower v sin / relative to high-density, intense 
star-forming regions. However, Huang and collaborators suggest that the lower v sin / 
in the field is simply linked to a higher average age for field star s, which may be dis- 
persed from cluste rs, reflecting an evolutionary spin-down effect ( Huang & GieslllOOSi : 
Huang et"al]|2010h . They demonstrate that both field and cluster stars show an anticor- 
relation between v sin / and log g, confirming this spin-down. 
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Figure 1. Preliminary distribution of vsin/ for solar neighborhood B stars 
(Braganga et al., in preparation). 

However, the field vsin / distribution strongly appears to be bimodal. We ob- 
served a large sample of solar neighborhood B stars with the MIKE echelle spectro- 
graph at Magellan, and removed known cluster members, yielding 240 remaining stars 
(Bragan9a et al. 2011, in preparation). Interestingly, this only heightens the bimodality 
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apparent in th e v sin / distributi on (Figure [p . The bimodality is also evident in the field 
star sample of ' Abt et all (|2002|) . selected by 'S trom et"aLl (l2005h and plotted by Wolff et 
al. (2007). There are also hundreds of stars in this sample, so the double-peak form 
of the vsin/ distribution is statistically significant, and similar to that in Figure [T] Gar- 
many et al. (2010) also identified a bimodal distribution in their sample of 150 Milky 
Way outer disk B stars, which were also observed with MIKE at Magellan. In all cases, 
we see a component of slow rotators which peak below vsin/ ~ 25 km s"^, and fast 
rotators with a much broader v sin / component peaking near 100 km s~'. Meanwhile, 
a sample of high-latitude, field main-sequence B stars studied by iMartinI (12006) does 
not show the bimodal distribution and apparently shows only the broad component of 
fast rotators. Martin argues that since these stars are outside the Galactic plane, they 
must be strongly dominated by runaway stars. Therefore, we suggest that the fast rota- 
tors seen in the bimodal distributions are runaway stars, while the population of slow 
rotators are field stars that formed in situ. This is consistent with the model of Wolff 
et al. (2007) that v sin / correlates with star formation density. Therefore, this may be 
strong evidence that massive stars do form in situ in the field. 



6. Runaway Stars 

Runaway stars ejected from clusters are commonly observed and understood to be a 
significant component of the field massive star population. O stars show a higher fre- 
quency, MwmdlO -30% runaways, relative to B stars, around 5 - 10% runaways 



(iMoffatll 19981 : lGiesl[l987 ). This mass dependence is an important clue to the origin of 



runaway stars, for which there are two general models. The first is dynamical ejection 
from 3- and 4-body interactions, with the latter taking place as binary-binary interac- 
tions, in which the ejected stars may or may not be separated. There has been substantial 
work done in simulatin g the se dynamical interac t ions, pionee red by, e.g., Poveda et aT] 
(ll967h . lMikkoial(ll983h . and Leonard & Duncanl(IT988Lll990h . The second mechanism 



to generate runaway stars is via supernova ejection, in which a binary companion ex- 
plodes and releases the surviving star as a runaway (e.g.. Zwicky 1957; Blaauw 1961). 

It is difficult to determine which of these mechanisms dominates, but that result 
would itself be a vital diagnostic for cluster dynamical evolution, binary fractions, and 
star formation. The strongest test is the frequency of runaway, double-lined spectro- 
scopic binaries (SB2) or visual binaries, which must be generated by the dynamical 
ejection scenario. Note that single-lined spectroscopic binaries may have a compact 
secondary, which could be generated by the supernova ejection model. Out of 42 run- 
aways in the Galactic field O star sample. Mason et al. (2009) identify 3 as SB 2 and 
11 showing visual multiplicity. The two categories are not mutually exclusive, and 
the visual multiples may include line-of-sight interlopers. Thus, the frequency of non- 
compact, multiple runaways may be around 10 - 20% of O star runaways. If this 
relatively high value is confirmed, it demonstrates the significance of the dynamical 
ejection mechanism. This may even dominate the massive runaway population, since 
this process also con tributes many single runaway stars. 



iBlaauwl (119931) suggested that stars ejected by the supernova binary mechanism 
should be He-enriched from contamination by the explosi on event. Such an enhance- 
ment indeed has been observed by Blaauw and by Hoogerw erf et all (1200 ih in a sample 
of local runaway stars at distances within 700 pc. Another diagnostic, proposed by 
Martin (2006), is that the runaway star's v sin / should correlate with its space velocity. 
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since the star's runaway velocity should be closely linked to its former binary orbital 
velocity. To produce a high, runaway velocity, the orbital velocity also must have been 
high, implying a tight binary that was most likely orbitally locked with the compan- 
ion, and hence the ejection and orbital velocities should be similar. However, Martin's 
(2006) sample of high-latitude field B stars does not show such a correlation, suggest- 
ing that the supernova ejection mechanism may not be important for that population. 
While Martin (2006) comes to the opposite conclusion from Hoogerwerf et al. (2001), 
he notes that the two studies are based on very different samples, with the latter per- 
haps representing a younger runaway population, owing to the stars' proximity within 
the Galactic disk. Further studies can examine the age-dependence of runaway frequen- 
cies and properties by targeting evolved stars (e.g.. lBerger & Giesll200ll : [Slettebak et al. 

Our understanding is still extremely uncertain, but the modest existing evidence 
described above perhaps slightly favors dynamical ejection over supernova ejection as 
the dominant mechanism for producing massive runaway stars. 
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Preliminary distribution of RIOTS4 radial velocities for SMC field OB 



Much more study is needed, especially observations to better quantify the proper- 
ties of runaway stars. With our RI0TS4 survey, we are able to measure radial velocities 
of our complete OB star sample with typical accuracy of 10 - 15 km s~'. Our prelimi- 
nary distribution (Figure|2ll has a FWHM ~ 60 km s^' , which includes large-scale SMC 
motions. In particular, the SMC Wing region is offset from the SMC Bar by about 30 
km s"^ From the fairly smooth velocity distribution, it is therefore difficult to infer the 
existence of any significant features that may originate from different runaway mech- 
anisms, but quantitative predictions are needed to fully exploit this complete da ta set. 
Our re sults are qualitatively similar to the distributions presented by Evans & HowarthI 
(l2008h . who compiled the radial velocities for samples of OB stars that include both 
cluster and field stars. Our sample is smaller than theirs, but with much stronger selec- 
tion criteria and spectral resolution that is almost double theirs. 



7. In-Situ Field Stars 



The existence of massive field stars that formed in situ may depend on the exact defini- 
tion of field stars (§1). While candidates for such stars do exist (e.g.. Seller et al., these 
proceedings), it remains to be determined whether any of these form in true isolation. 
But in any case, the sparsest regime sets strong constraints on the physics of star for- 
mation. For example, it is much m ore difficult to gen erate isolated massive stars with 
comp etitive accretion models (e. g., lBonnell et aLll2004i) than with core collapse models 
(e.g. jKrumholz & McKe"3l2008h . 
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There do appear to be examples of massive stars in sparse groups containing only a 
few other low-mass stars. Such groups might be remnants of rapidly-dissolved clusters, 
but it remains like l y that these are examples of extreme, low-density star formation. 



de Wit et all (|2004 120051) find that locally, about 12% of Galactic O stars are in such 



sparse groups. Using HST/ ACS observations, we found three SMC field O stars to be 
withi n sparse, low-mass groups, out of six apparently-isolated field O stars (Lamb et ahl 
We carried out Monte Ca rlo simulations based on the -2 power law clustering 
(equation [T]) and a ' Kroupal ([2001') IMF, assuming these two relations to be completely 
independent; in particular, we assumed no dependence of the stellar upper-mass limit on 
cluster mass. The observed sparse O star groups fall easily within the allowed parameter 
space of these simulations when examining the mass ratio of the two highest-mass stars 
and the maximum-mass star m^ax per cluster. However, these sparse O star groups are 
inconsistent with the proposed relation between mmax and cluster mass, suggesting that 
in general, the integrated galaxy IMF (IGIMF) should not be steeper than the IMF in 
individual clusters. 



8. Conclusion: In-Situ Field Stars In-Situ and Runaways 

The evidence presented above supports the existence of both runaways and an in-situ 
component contributing significantly to the field massive star population. We com- 
monly observe runaway stars and we know that they exist; thus, the main question is 
whether field massive stars also form in situ. I suggest that several lines of evidence 
support their existence. First, there are many instances of observed high-mass stars that 
are in sparse, low-mass groups, or even candidate in situ isolated stars. Second, the 
bimodal v sin / distribution strongly suggests that the slow rotators are in situ field OB 
stars. Third, the upper IMF in the field appears to be significantly steeper than mea- 
sured in clusters; whereas, if the field consists purely of runaways, then its upper IMF 
should be flatter than in clusters, since O stars have a higher runaway frequency than 
B stars. And fourth, there is no discontinuity between the field stars and clusters in the 
cluster mass function (equation [T|). 

The runaways originate from either dynamical ejection or supernova ejection, and 
presently, the evidence does not strongly favor one scenario over the other. There may 
be slightly stronger support for dynamical ejection, since we do see a non-negligible 
frequency of double-lined spectroscopic binaries and visual binaries, both of which 
can only originate by dynamical ejection. The lack of correlation thus far, between 
rotation and space velocity also suggests that supernova ejection is not a dominant 
process. However, observed He abundance enhancements in runaway stars do support 
the supernova ejection model. All these results are based on limited, individual studies 
and far more work is needed to establish the trends, and better probe diagnostics of the 
ejection mechanisms. 

Our understanding of the field massive star population is gradually emerging, but 
much more study is needed to fully illuminate this enigmatic, but substantial, popula- 
tion. An especially critical area is the binary population: some of the most powerful 
diagnostics will be revealed by statistics of the binary orbital parameters. These are dif- 
ficult and time-consuming to obtain, but they will be powerful diagnostics of runaway 
ejection mechanisms, constrain star formation and IMF theories, and cluster evolution, 
as well as the nature of the field population. Binary population synthesis models and N- 
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body dynamical simulations are already mature areas of study, and more observations 
in particular are needed to test the models. 
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